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SOUND RADIATION FROM A HIGH-SPEED AXIAL-FUOW FAN DUE 
TO THE INLET TURBULENCE QUADRUPOLE INTERACTION 
by Marvin E. Goldstein, BurtM. Rosenbaum, and Lynn U. Albers 

Lewis Research Center 

SUMMARY 

One of the dominant sources of broadband noise (and possibly also of pure tones) in 
aircraft engine fans is the turbulence drawn into the fan inlet. A number of investiga- 
tions have analyzed this noise source by using a dipole model. However, at the high 
subsonic Mach numbers at which some current engine fans operate, the quadrupole 
mechanism may be dominant over the dipole mechanism. In such a case, analyses based 
on the dipole model are unsatisfactory. No analysis of this noise source which is based 
on the quadrupole model and is suitable for high subsonic Mach numbers has been given. 

The present report is written to fill this gap. 

i ■ 

A formula is obtained for the total acoustic power spectra radiated out the front of 
the fan as a function of frequency. The formula involves the design parameters of the 
fan as well as the statistical properties of the incident turbulence. Numerical results 
are calculated for values of the parameters in the range of interest for quiet fans tested 
at the Lewis Research Center. As in the dipole analysis, when the turbulence correla- 
tion lengths become equal to the interblade spacing, the predicted spectra exhibit peaks 
around the blade passing frequency and its harmonics. There has recently been con- 
siderable conjecture about whether the stretching of turbulent eddies as they enter a 
stationary fan could result in the inlet turbulence being the dominant source of pure 
tones from nontranslating fans (i. e. , fans in nonmoving aircraft). The results of the 
current analysis show that, unless the turbulent eddies become quite elongated, this 
noise source contributes predominantly to the broadband spectrum. 


INTRODUCTION 

The noise generated by single-stage subsonic fans consists of a number of discrete 
tones (predominantly at the blade passing frequency) superimposed on a broadband spec- 
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trum. It is generally accepted that the broadband noise is generated by such factors as 
interaction of the blades with patches of inlet turbulence, irregular vortex shedding by 
the fan blades, and random modulation of rotor wakes impinging on a downstream stator. 
The discrete tones are believed to be caused by periodic spatial nonuniformities entering 
the fan and by the passage of its wakes across the stator. However, if the blade passing 
frequency is large compared with the characteristic turbulence frequency, the sound 
generated by the inlet turbulence interaction will be concentrated around the harmonics 
of the blade passing frequency and the spectrum will appear to contain tones of finite 
width. Thus, the inlet turbulence interaction may also contribute to the pure tones. 

This report is concerned with the sound generated by the interaction of inflow tur- 
bulence with a high-speed compressor fan. The turbulence in the flow entering a com- 
pressor fan (both in the casing boundary layer and in the main stream) can produce noise 
through two mechanisms. One is a quadrupole-like mechanism which is the result of the 
fluctuating Reynolds stresses in the volume exterior to the blades produced by the inter- 
action of the turbulence with the potential -flow field of the rotor. The other is a dipole - 
like interaction which is the result of the fluctuating lift on the fan blade resulting from 
the fluctuations in angle of attack caused by the turbulence. 

Until fairly recently, all analyses of fan noise considered only the dipole mechan- 
ism. However, at the operating Mach numbers of modern fans the quadrupole interac- 
tion is probably dominant. 

Analyses of the inlet turbulence noise generated through the dipole mechanism have 
been given by Sevik (ref. 1) and Mani (ref. 2). The possibility that this noise is gener- 
ated through a quadrupole mechanism was first pointed out by Ffowcs -Williams and 
Hawkings (ref. 3) and was later investigated in detail for a low -speed fan by Chandra- 
shekara (ref. 4). Chandrashekara carried out a combined theoretical and experimental 
investigation of a low -speed (tip Mach number less than 0. 3) unducted fan. He compared 
the "tone" noise predicted by both his analysis and the dipole analysis of Mani (ref. 2) 
with his measurements and found that the level of the tones predicted by the quadrupole 
model was considerably below the experimental values, while the dipole model predicted 
the correct level. However, the ratio of the quadrupole noise to the dipole noise should 
vary as Mach number squared and should increase with increased blade loading (ref. 4). 
Hence, at the high tip speeds and blade loadings at which current jet engine fans oper- 
ate, it is likely that the quadrupole mechanism will dominate. However, even though 
inlet turbulence is a dominant noise source in fans, no quadrupole model suitable for 
high subsonic Mach numbers is available for calculating this noise. The present report 
is written to fill this gap. 

In order to construct a model which applies to a ducted fan at these high speeds, it 
is necessary to use a model which differs from the one used by Chandrashekara in a 
number of important respects. First, the present model considers a fan operating in an 
infinite duct containing a mean axial flow and not in free space with zero mean flow. In 
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Chandrashekara's analysis the variation in retarded time across the turbulent eddy was 
neglected. However, as pointed out by Ffowcs-Williams and Hawkings (ref. 3), this is 
not a good approximation at high tip Mach numbers, especially for the 11 tone" noise, 
since the sound is concentrated at the blade passing frequency, which is quite high, and 
hence corresponds to a condition where the wavelength of the sound can be smaller than 
the "eddy” size. Thus, this omission also serves to limit Chandrashekara's analysis 
to low Mach numbers. The present analysis also explicitly accounts for the convection 
of the eddies by the mean flow, which from studies of jet noise is found to be important 
for properly modeling the turbulence. In the present analysis the effects of compres- 
sibility on the rotor potential -flow field are included. This is important because of the 
high Mach numbers of interest. A more detailed model of this flow field is also adopted. 
In real turbulent flows the turbulence is not isotropic, but the eddies tend to be elon- 
gated in the direction of flow. Since the flow into a stationary inlet involves a large con- 
vergence of the streamlines, it has been argued that the turbulent eddies can become 
quite elongated in the direction of flow. We have therefore adopted a turbulence model 
which is general enough to account for this effect. 

The analysis leads to a formula for the power spectral density of the sound field 
emitted through the fan inlet as a function of various aerodynamic and geometric param- 
eters of the fan and the properties of the inlet, turbulence. This formula cannot only be 
used for incorporating the quadrupole inlet turbulence noise into fan noise calculation 
packages and acoustic liner design programs, but should also be useful for making pa- 
rametric studies of the sound field produced by this noise source. Hence, numerical 
results are computed for parameters covering the design range of the various Quiet Fans 
which have been tested at Lewis. The predicted spectra exhibit a number of features 
which are observed in the spectra of actual fans. But it is found that the width of the 
tones predicted by the analysis will generally be broader than those observed in the fan 
spectrum unless the turbulent eddies are quite elongated (typically of the order of a 
couple of feet). However, there is mounting evidence that the turbulent eddies drawn 
into a nontranslating fan do indeed become very elongated and as a result are the domi- 
nant source of pure tones in turbojet engines on the ground. This elongation does not 
occur during forward flight and, in. fact, recent measurements by Cumpsty and Lowery 
(ref. 5) indicate that the tones produced by aircraft engine fans are considerably reduced 
in flight. The present analysis should be helpful in resolving this question. 

ANALYSIS 

Basic Formulas for Quadrupole Sound Emission 

It is shown in reference 6 that the density fluctuation p'(x, t) at the point x and the 
time t due to the quadrupole radiation from a fan in an infinite duct containing a uniform 
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axial flow with velocity U is governed by the equation 



9G Tjjdydr 


3y_j 

where Cq is the speed of sound, T is the time, 

T^(y, T ) = pvjvj + 6.j(p - CqP^ - 
is Lighthili’s stress tensor based on the relative velocity 


u) 


e ij 


( 2 ) 


V 1 5 v i * 6 li U 


( 3 ) 


instead of the actual velocity v^, p is the pressure measured above the ambient pres- 
sure p 0 , and is the viscous stress tensor. (All symbols are defined in appendix B. ) 
The volume p{r) over which the integration is carried out is the entire region of the duct 
external to the blades. The uniformly moving medium, outgoing -wave, infinite -duct 
Green’s function G is given by (ref. 6) 


G(y, t x, t) = 


±Z 


V r/ x 2’ x 3> 


m, n 


m, n 



where k = u>/c Q , M = U/cq, j 8 = \ 1 - and 


4> m n denotes a doubly infinite set of 
eigenfunctions with eigenvalues tt m n determined by solving the Helmholtz equation 


+ k & - 0 

m, n > m, n m, n 



( 5 ) 
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on the cross-sectional area A of the duct subject to appropriate homogeneous boundary 
conditions on the duct surface. Thus, when the surface of the duct is rigid, we must 
require that 

3 $ 

— = 0 on S (6) 

3n 

where 3/3n denotes the normal derivative to the duct surface S. The remaining quan- 
tities in equation (4) are 


r m,n s / I*m,nl 2dy 2 d y3 (7) 

and 

where, in order to ensure that G represents the outgoing-wave Green T s function, we 
must choose the branch out of the square root to be as shown in figure 1. 

When the observation point x is far enough from the source (i. e, , in the far field), 
the pressure and density fluctuations will be small and we can use the linear approxima- 
tion 


Rectangular Duct Model 

In order to model a real fan, the most appropriate cross-sectional shape of the duct 
is an annulus such as that shown in figure 2(a). The eigenfunctions 4> m will then be 
combinations of Bessel functions. The analysis can be considerably simplified, how- 
ever, by assuming that the duct is "unrolled" into the rectangular strip shown in fig- 
ure 2(b). (The larger the number of blades, the more closely the rotor approximates a 
two-dimensional disturbance pattern with subsonic phase speed. ) The width 6 of the 
rectangular duct is equal to 2n times the mean radius R of the annular duct, and its 
height b is equal to the outer radius minus the inner radius of the annular duct. In 
this case we must require that 4> m and its normal derivative take on the same values 
on the surface S a as they do on the corresponding points on the surface (periodicity * 
conditions). It is easy to verify that the eigenfunctions of equation (5) which satisfy this 
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condition together with the boundary condition (6) on the upper and lower walls are given 
by 


where the y 
ar e given by 


Hence, 


and 


G(y, t x, t) = 


i27rmy 9 /6 nay« 

= e 6 cos for 

m,n b 


n = 0, 1, 2, . . . 


,ni - 0, ±1, ±2, 


•coordinate system is shown in figure 3 and the corresponding eigenvalues 


K m, n 



( 9 ) 



(10) 


i 

2x6b 




n7iy 3 

cos 

b 


mrxo 

cos 

b 


1 + 6 


n, 0 


i2ma(y2-X2)/6 
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where the upper sign corresponds to upstream propagation (x^ < y^) and the lower sign 
corresponds to downstream propagation (x^ >yi>- 


Transformation to Rotor Coordinates 

In order to make the limits of integration in equation (1) independent of time, we 
introduce the coordinate system by 


>?i = Vi + 6 i2 U t T (12) 

where U t denotes the velocity of the blade row and is taken as positive in the direction 
shown in figure 3. But since for each value of (y^, y^) T?j takes on the same value at 
yg = -6/2 as it does at y2 « +6/2, it can be extended to a continuous periodic function 
from to +°°. And since the Green f s function is also periodic, the integration with 
respect to can be carried out over any region vq of fixed shape which coincides 
with v ( t ) at some definite instant of time. Hence, equation (1) becomes 


p(x, t) = 




8r ?i 


T)dq d T 


where the linear approximation (8) has been used* 


(13) 


Derivation of Spectral Equations 

Rather than deal with the pressure directly, it is convenient to deal with its spec- 
trum P(w) defined by. the generalized Fourier transform 



(14) 


It then follows upon inserting equation (12) into equation (11) and using the result in equa- 
tion (13) that 
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(15) 


oo 


P( “ ) = £ Z P m ; n^ c m,n® 


m--°° n=0 


where 


, - f 

^0 


,2 i 


„ - Wh “ w - w 


/ 


i(w-2m7rU f /6)r 

i L 'T'T 


T-. dr dr; 


(16) 


ff rV Mk± V^ 


t*2 = 2m7r/6 


ffg = n?r/b 


(17) 


cos I 


MlTTX, 


b 


1 V 2m7rx 9 


m, n 


2^6b (1 + 8 ^ 


exp<-i 




(18) 


and the plus sign corresponds to upstream propagation while the minus sign corresponds 
to downstream propagation. 


Radiated Acoustic Power 

The quantity which is perhaps of most interest is the total acoustic power per unit 
frequency radiated (upstream/downstr earn). This can be calculated by integrating 
the axial component I^x) of the intensity spectrum over the cross-sectional area of the 
duct. Thus, 1 


*The factor 2 arises because both 1^ and I contribute to the acoustic power 
with frequency co. Since 1^ arises as a Fourier transform, it is defined for positive 
and negative frequencies. But since 1^ = I and the experimental power is defined 
only for positive frequencies, the factor of 2 must be included. 
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(19) 


^(w) 


/ '6/2 /-b _ _ 

lim 2 / / I, ,(x )dx„ dx„ 

x r » 7-5/2 7 0 <*> 2 3 


We shall suppose that the sound field is time stationary. Then it is shown in reference 6 
that 


I ^ lim 
w X— 00 T 


(l + m 2 )pv + -M- | p | 2 + 

p 0 c o 


p 0 c 0 M |v 


( 20 ) 


where T is the integration time for the Fourier transform and V is the generalized 
Fourier transform of the axial acoustic velocity v|(x,t). Thus, since p and vj obey 
the linearized acoustic equations at large distances from the fan 


Hence, 


3p =-Po c of ik - M r-V 


dx 




dx. 


Inserting equations (15) and (17) into this relation shows that 


QO CO 


V = — T y n P m ,»C* „(*) 

i—t m,n m,n' ' m,n v ' 
p 0 c 0 m=-« n=0 


(21) 


where 


Mk ± 


n 


m, n 


k±Mk m,n 


( 22 ) 


Since equation (18) shows that 


■6/2 rb _ 6 

_C„ „ dx 0 dx„ = BiRJhlL 


ro/z r 

J-fr/9 


■6/2 -m) m, n p, q 2 3 , 2 

7 8jT 2 6bk^ (1 + 6 n ) 

m, n x n, U' 
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Substituting equations (15) and (21) into equation (20) and using the result in equation (19) 
show that 


^ s± (o>) 


T/3 4 k 



ft 


m, n 


‘"‘W K * S ". °>V «( k * “V 

all n==0 with 

k2> ^ 2 4,„) 


(23) 


where, since cutoff modes do not contribute to equation (23), the sum in this equation is 
carried out only over propagating modes and 


ft* „ s lim |p± ( w ) 

m, n T _ 0O ,p ' m, n v 7 


Inserting equation (16) into this relation shows upon using the Wiener- Khinchin relations 
for stationary processes 


ft! _ = 2 tt 
m, n 



exp ji ~ ’^1^ + a 2^2 ~ ^2^]} cos ^ a 3^3^ cos ^ Q, 3 r? 3^ 
\ d ?l' k dv l 


where 



-i(cu-*2ni7rU t /a)T 

dr dr) &}' 


(24) 


<* m (?, ?, t) s 4 .® t)T^ (?, t + r) 


(25) 


and the overbar denotes the time average 


f s lim 
T 


ml/ 

«fcOO l T' — T 


T/2 

T/2 


f(t)dt 


(26) 
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Source Model 


In order to use these results to predict the sound emission, it is necessary to deter- 
mine the source term ^ , As in any aeroacoustic problem, this is extremely diffi- 
cult to do exactly since it requires a knowledge of the complete flow including the acous- 
tic field. It is therefore customary in aeroacoustic problems to develop an approximate 
model for the source term. In this section the general structure of this source model 
will be evolved, and it will be used at the end of the section to simplify the equations de- 
rived in the preceding section. The source model will involve contributions from the 
turbulence and the phase -locked rotor velocity field. In the next section a model for the 
rotor velocity field will be developed, and the results will be used in subsequent sections 
to simplify the formulas for the sound field. A model for the turbulence correlation ten- 
sor will then be introduced, and with this in hand the final equations will be derived in 
the last section. 

As a first step in evolving the source model, we shall replace Lighthill's stress 
tensor by the Reynolds stress (ref. 7). (This procedure is followed in most anal- 
yses which do not involve combustion. ) 


T Jj = po v M (27) 

We shall further assume that the relative velocity vj can be decomposed into the sum of 
a random fluctuating velocity u^ (with zero mean) associated with the inlet turbulence 
and a velocity w i associated with the phase -locked rotor pressure field. Then w i is 
independent of time in the if- coordinate system and 

vt(v, *) = u i(v, t) + Wj(^) (28) 


= 0 (29) 

where u. and w- denote velocities relative to a coordinate system moving uniformly in 
the y^ -direction even though they are expressed as functions of the rotor -locked 
rf-coordinate system. 

Inserting equation (28) into equation (27) and using the result in equation (25) shows 

that 
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+ R ij, k w l + R ij, l w k + R i kZ w j + R j Zk w i + Time ' inde P end ent terms 


(30) 


where 


R ijkZ = u i u j u W 


R ij, k = u i u j u k 


R i,ik = Vj u k 


are turbulence velocity correlations and the unprimed quantities are evaluated at r\ 
and t while the primed quantities are evaluated at rp and t + r. 

When the integral with respect to T in equation (24) is carried out over the time- 
independent terms, it wi ll yield a term proportional to 5(w - 2m7rU t /5). But as long as 

relative Mach number + (U t /c 0 ) 2 is less than unity, k 2 will be less than (3 2 K m 

with cu = 2mffUj./5 and all modes will be cut off. Hence, for subsonic rotors the time- 
independent terms will make no contribution to the sound field and can be omitted. 

We shall assume (as is usual in compressor analyses) that the phase-locked rotor 
flow field is two dimensional so that Wj is independent of y^ = ?jg. It is convenient to 
introduce the new variables 


I s 


Vi 1 




( 31 ) 


We shall no w^ suppose that the turbulence is homogeneous. Then (considered 

as a function of £ and Y) will be independent of Y^. We shall also suppose that the 
region of integration can be extended to the interior of the blade row and that the turbu- 


lence correlation length in the y 2 ~ and y 3 -directions is small compared with the lengths 
of <5 and b in those directions- Then the integration with respect to If can be carried 
out over all space and the integration with respect to Y carried out over i/q. (Hence, 
upon using the identity) 


cosf- 


rur \ -_„/n7r » \ 1 


J njr|o 2nj7Yo N 

'cos + cos 

lb b 


equation (24) becomes 


ft x 
m, n 



d T d£ dY 


(32) 


where the notation / indicates that the integration is to be carried out over all space. 
Now the only term depending on Y g is cos(2n7rY 3 /b), and the integration with respect to 
this variable is from 0 to b. But since cos(2njrY^/b) and all its derivatives integrate to 
zero over this range whenever n * 0, we can carry out the integration with respect to 
Y 3 in equation (32) to obtain 


ft 


m. 



4 -i (of| i+af^) 


a’e ' * '* “ cos 


L 


-i(co-2m7xU,/6)r ^ 

X I e ^ ijkz d T d£ dr^ (33) 


where Aq denotes the cross-sectional area of any plane, y 3 » Constant, included be- 
tween the duct walls and the exterior of the blades. 
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The first term on the right side of equation (30) represents the effect of the interac- 
tion of the turbulence with itself. This noise generation mechanism is the same as that 
which produces jet noise and is known to be a highly inefficient producer of sound rela- 
tive to the other noise mechanisms which are present (refs. 3 and 7). We shall there- 
fore neglect this term. When equation (30) is substituted into equation (33), we therefore 
obtain after interchanging dummy indices 


fi m,n = 4 ^ 1 + 6 „,0> 


4 -i a'iq+Q' 2 ± 

3 e cos q-oIq 

3«i a« k 


X 


*/ 


R ik e 


-i(o)-2m7rU 


,/«7 , r 

a ”i"i / "i, 

«/_oo 


-i(w-2mjrU./6)r 


kZ 


l/-0O 


-i(oJ-2m7rU t /6)r 


di 


dr 


where the rotor velocity integrals H^, H-, and H! are defined by 


H 


\L s [ w j w l A *2 ir >l 


H j a jf W j ^2 dr *l 

A o 

H r/ 




dr 


(34) 


(35) 


(36) 


Rotor Velocity Field Model 

We shall now develop an approximate model for the phase-locked rotor velocity field 
Wj. Most fan blades have fairly small camber, and linearized theory can be used. But 
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since the quadrupole source is of interest at high subsonic Mach numbers, we must take 
into account the effects of compressibility. It is convenient to introduce the X- 
coordinate system, which is alined with the blades as shown in figure 4. It is well known 
that the potential -flow field about such a blade row can be represented by a uniform ve- 
locity plus a distribution of line vortices along the blade chords. In order to simplify 
the analysis, we replace this vortex distribution by a discrete set of concentrated vor- 
tices lying at the midchord point. Let V-^ and Vg denote the velocity components along 
the X^- and Xg -directions, respectively {relative to the blades). Then according to 
linearized theory the velocity induced by a single compressible vortex of strength Tq at 
the origin is given by 


where 


and 




X 1 + ^r X 2 


(38) 


Hence, the concentrated vorticity approximation to the velocity field about the cas- 
cade is given by 



(39) 


where Vq is the uniform (in general, complex) velocity (which will be chosen to make 
the flow far upstream of the cascade equal to U ) and 

A 0 = A(/3 r cos y - i sin y) (40) 


where A is the interblade spacing and y is called the stagger angle. But upon using 
the relation (ref. 8) 
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oo 



Z - inn - 


= coth Z 


equation (39) becomes 


But since 


we must put 


Hence, 




— coth 


2A 0 i 


lim coth -1 

X l--°° A q 


v 0 


r 


o 


^r 


2A 


0 


V, 


Vi-i 


= u r + 


1 + coth 


ttZ 


i2A, 


(41) 


It can be seen from figure 4 that the X^, coordinate system is related to the 
rotating duct alined coordinates by 


X 1 " ^1 cos y + ? ?2 sin y 


X g = - 7 ]^ sin 7 +t ?2 cos y 

Hence, it follows from equations (38) and (40) that 
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z 7 ?1 (cos y - i p r siny) r/ 2 

Aq A(j3 r cos y - i sin y) A 


*1 


/3 r + i 


(l - j3 2 )sin y cos y] r) 2 


A^ cos 2 y + sin 2 y 


) 


But upon using the velocity triangle in figure 4, this can be written as 



(42) 


where 


z = ?)-|D + i?72 


(43) 


\ 


P + iMM. 



(44) 


and 



(45) 


is the Mach number of the blade row. 

The velocity components w^ and w 2 are in the tj^- and tj 2 -directions and are 
measured relative to a reference frame moving with a uniform velocity U in the re- 
direction, while the velocities V x and V 2 are in the X^- and Xg-directions and are 
measured in a reference frame moving with the cascade velocity U t - Hence, it follows 
from figure 4 that 

V, = (w x + U)cos y + (w 2 + U t )sin y 
V 2 = (w 2 + U t )cos y ~ (w x + U)sin y 

= w 2 cos y - sin y 
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Then 


/ V, \ - iw„D 

(0 r cos y - i sin yHVj - i— - U r j - -i — 

\ / ^r 


Multiplying equation (41) through by (j8 cos y - i sin y)0 r and using equations (40) 
and (42) show that 


p Q 

j3 2 (w 1 - iWpD) « -5-1 A + coth 

2 2iA y A y 


(46) 


Then 


r 

lim iS 2 (w 1 - iw 9 D) = — — 
7?l ^ + oo z 2iA 


Equating the imaginary parts of this expression shows that Aw 2 , the change in w 2 
across the blade row, is related to Tq by 



where 9 = -Awg/Uj. is called the work coefficient of the ran. Hence, equation (46) can 
be written as 


9 _ iOB U, / \ 

P&l - iw 2 D) - — U + coth — j 

i^ r U t e 772 / ^ 

2 sinh ~ 

A 


By using the geometric series 


1 

1 - z 



Z P 


for z | < 1 
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(with z — exp(±2ffz/A)), this becomes 


^ 2 (w 1 - iw 2 D) = 


Equating real and imaginary parts now shows that 

w 1 i3 2 - w 2 M t M ^ ( e 2 ^P z /A _ e 2j7pz/Aj f 0r ^ < 0 

2 p=l * 

w 2 = £ (e 2irpz/A + e 2irpz/A ) for ^ < 0 (47) 

2 P=1 

Or adding the results 

(De^^-De 2 ^^) for ^ < 0 (48) 

2 P=1 

The noise generated by inlet turbulence probably has its largest effect on the sound field 
passing through the fan inlet. We shall therefore restrict our attention to the upstream- 
propagating waves. We shall also suppose that only the front half of the cascade 
potential -flow field is effective in interacting with the inlet turbulence and generating 
sound which propagates upstream. Hence we set 

W 1 ~ w 2 ~ 0 for ^>0 (49) 

Evaluation of Rotor Velocity Field Integrals 

The results of the previous section will now be used to evaluate the integrals defined 
by equations (35) and (36). It follows from equations (47) and (48) and the orthogonality 
properties of the complex exponentials that 


~ iU t'V]C e 27rpz / A if ?7 X < 0 

p=l 


il W E e_2,ipz/A if ^>0 

P=0 
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(50) 


Hj = Hj = 0 


Next (in order to facilitate the analysis) we define the rotor velocity correlations by 

. 6/2 

■ 6/2 


.6/2 

LjZ s / w,(5)w i (^+ 


(51) 


Equation (35) shows that 


H 




/oo 

L ji dJ )l 


(52) 


provided we extend th^. region of integration to the interior of the blades. Inserting 
equations (47) and (48) into equation (51) now shows that for 

il j < minHpO} 


L 21 = 


6U?f> 2 


re \ ' „ /2 b \ 

i— At > Dexp 2jEf-I£L +D « 1+ i{. 

2 A Vs 2 / 

P =1 L p 


'12 


6U?0 2 

-^-AtB \ 

p=i 


exp 


2np 


— W + D h + i h 


L 22 " 


2 2 00 
6U“ir 


ft V exp 

4 A lV' 


Vi + D ?1 + i?2 


and 


. , 2 1-Mf 
L U -| D | L 2 2 T~ L 22 

F 


Lji =o 
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otherwise. Inserting these results into equation (52) now shows that 


1 - M? 

h u = — h 22 



Development of Equation for Upstream Sound Propagation 


Inserting equation (50) and equations (53) to (55) into equation (34) (with the plus sign 
only) now shows that 


a* =y o' 

m, n m, n 

P =i 


(56) 


where 


„ 2rrp?6bU?fl 2 

a? 


m. n 


0, 


(i 


o / £ i + a 2^2^ i r "ib'qlo i(ttolo+6ojJT+6ob.^)l 

+ ^.,0>Wk J e 1 1 2 2 i|e 3 3 + e 3 3 31 3k J 


/ 


X I R ik e 


-i(ai-2mjrU t /6)r A -((3 //r)U 2 |(2irp/A) 

e cos 

47Tp 




Lne ! 


d£ dr 
(57) 
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and 


®0 = a \ + a 2 ' ( M «2 ' M t a l )2 (58) 

In order to develop a model for the turbulence correlation tensor R ik , it is conven- 
ient to introduce a coordinate system which ’’moves with the turbulent eddies. ” Thus, 
we introduce the £ -coordinate system by 


Ci = «i - V 

c 2 = - V 

C 3 = *3 


where U denotes the convection velocity of the turbulent eddies. It is not necessarily 
c 

equal to the axial velocity U. Equation (7) can now be written as 


ft 


2 7TPq tt2 * 2 


m, n 


U 2 i> 2 6b(l + 6^ o)«Q« 1 a , k r-Kaj^+OjC 2 ) l ^ 

xf ^ ex Pp<“ + "lV - + VI ^ 


-ia 3 ? 3 l(o 3 ?3+6 3i iT^3 k j)l 




fMM t 

7" 


(U c r + ?i) + ?2 + % 


kd T 


(59) 


where 


# ik (f, t) = R ik (|, r) 


(60) 


is the moving-axis turbulence correlation tensor. 


Model for Turbulence Correlation Tensor 

In order to proceed, it is necessary to determine the turbulence correlation tensor. 
We first assume that the moving -axis turbulence correlation tensor factors into space - 
and time-dependent parts. Thus, we put 
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( 61 ) 


<% = 
lk" 




where is the decay time of the moving -frame turbulence correlation, 

simplest assumption which can be made is that the turbulence is isotropic, 
general form of which satisfies the requirements of continuity (i. e. , 

matically possible) is 


Perhaps the 
The most 
which is kine- 


1K [ 2 d^y 1K 2£ 3£ 1 K 


where f can be any function of the magnitude £ of £ only. A reasonable assumption 
in the moving frame is that 


f = u 2 e -^ 2 


(62) 


where l is the turbulence correlation length. The turbulence in the inlet of a real fan 
is probably not at all isotropic, with the eddies probably being elongated in the direction 
of flow. In order to develop a kinematically possible turbulence model with scale an- 
isotropy, we follow Kraichnan (ref. 9) and note that the correlation tensor is a 

contravariant tensor of second rank that satisfies the continuity condition 

a 

= o (63) 

3 ?k 


which is a tensor equation. Then corresponding to the transformation of coordinates 


~ -1 

q = r h 




«a 


«2 


?3 


^3 


> 


(64) 


the transformation 
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-1^0 


Ik 


for k * 1 > 


ft - gp 

ik ^ik 


for i, k * 3 


J 


(65) 


ensures that if obeys the equation 3^?^/ 3?^ - 0, then ^9^. will obey equation (63) 

and hence be kinematically realizable. Thus, if we put 




V ] 


- c 2 /* 2 


( 66 ) 


the correlation tensor ^ k (C) will be kinematically realizable (i. e. , satisfies continu- 
ity) and will represent a turbulent flow where correlation length in the axial (y^) direc- 
tion is r l while the correlation length in the transverse direction is l . 


Rearrangement of Equation for Upstream Sound Propagation 
Substituting equations (61), (65), and (66) into equation (59) yields 


ft- 


p . 2TPjj6bufufo 2 (l + 


m, n 




x § e 


-i(o' 1 C 1 +a 2 C 2 ) _i «3?3 i t ff 3^3 +6 3i 7r+6 3k ,r ) 


-le 
2 


+ e 


i* 


(67) 


where 


s A_ 

m > n 4 np 


J exp p((D + oqtyr - w f | T | - + u < 


,r|^e 



r r 

- cos< 

See! 

A 


MM 


^(UcT + Ci) + ? 2+ U t r 


(68) 
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'rot- for i = 1 


L.ar. for i * l 


( 69 ) 


Equation (67) can be written as 


ft p M = lim 
m ’ n X— 1 


2^6bU 2 ufe 2 (l + 5 ni 0 )ag 
p r 2 


.iji + s a (i-.t-L'j 

j2 m2 \ 2 3^7 


^ ax 


J'^-i{a 1 i; 1+ a 2 Z 2 )X i / -i a 3 ? 3 X + ^ V ? 3 \ e -? 2 /2 2 M P 


where 


~2 2 2 2 2 
a - r dj + a 2 +. a 3 


(70) 


Upon carrying out the integration with respect to this becomes 

,2 EUTT 2..2 rt 2/i c \ 2 

r» 

ft 


_■ 2ip?6bU 2 u?6 2 (l + 6„ Jan 

jP = _2 Li Ss-5 — 2 ^jr lim 

m ’ n „ _2 A-l 


^r r * 


J - JL + « Vi - i 2 V 


Z 2 3A 2 


2 cZ 


Z e 


-(» 3 Ai/2)‘ 


a 
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x / / exp - i(a 2 ^ 2 + cy 1 ? 1 )X - 




M m,n d ?l d ?2 ‘ 71 » 


It is shown in appendix A that 


M p 
m, n 


tfip/AHSj , -(2irp/A)i5 2 ■ 

Vn,p^l ) + e in, n, p ^1^ 


where ICl are given by equation (A2). Inserting this into equation (71) and carry- 
P . 

ing out the integration over ^ shows that 
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r 2 pgcMflft a (i + a n< p) j2 


a 


m. n 


V' 


lim 
A— 1 


.iJI + tf/i-iJA 

L i 2 ax 2 V 2 5l { 


2 ”(®|j^/2) 

l e 


r ? 2 
e 


-2 2 (a 2 A-2rrp/A) 2 /4 f -«,/r 0 2 -la,?,X 

/ e IC _ _ d£< 

J m, n, p * 1 

-2 2 (ff ; ,X+2irp/A) 2 / 4 f . ~1 

+ e 2 y e Vn,p^lJ 


( 72 ) 


But 


/- 


-? 2 /r 2 2 2 -i<y 1 X? 1 


^n.p^l 



where c. c. denotes the complrex conjugate of the preceding expression. A similar ex- 
pression holds for the integral over K It can be obtained from the preceding 

m, n, p 

equation by replacing p by -p. Hence, equation (72) can be written as 
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ir 2 pg6bU 2 u 2 fS 2 ^(l + 6 ) , 

: — a 0 

c 0^r r 


Q p 
m, n 


x lim 
A— 1 


iii + S 2 fi 

l 2 d A 2 ' 


1 j\ 

2 az/ 


o -(a<,AZ/2)‘ 
l 6 e 6 


$ e Jr p + r” p ] 

Ke \ m, n m, nj 


where 



where is the same as defined by equation (17), W(Z) is defined in refer- 
ence 9 by 


and 


W(Z) a e Z erfc(iZ) 


sgn(x) = 



if x > 0 
if x < 0 


Derivation of Final Sound Power Equation 

Substituting equation (73) into equation (23) (with the upper sign) and using equa 
tions (17) and (58) show that 



(15) 


P 0 ugufo 2 V^ 

4c 2 /3 r r 


^ + (co) = 


z 


(All m and 
n > 0 with 


1 + < n - ( M - M tV / 

* t m, n 


k 2 > 


n) 


X lim 
A— 1 


1 



z 3 e -(nj7AZ/2b) 2 



m, n 


where 


m, n 


is given by equation (74) 


2mnp 



n 



(76) 


and equations (17) and (70) show that 


a 


2 


(Mk + n ) + (■ 

P V 


2m n 


6 / 


\ +/52 

b 


( 77 ) 


DISCUSSION OF RESULTS 

Equations (74) to (77) are the final results. They can be used to calculate the power 
spectral density of the quadrupole sound emitted by a fan as a result of inlet turbulence. 
In order to do this, it is necessary to specify the longitudinal correlation length l of 
the turbulence in the direction transverse to the flow, the ratio r of the longitudinal 
correlation length in the axial direction to the transverse longitudinal correlation length, 

the mean square turbulence velocity u^, the turbulence convection velocity U , and the 
characteristic frequency of the moving-frame correlation. The parameters b, B, 

A, and 6 = BA characterize the geometry of the fan. The remaining parameters depend 
on the operating conditions of the fan. They are the tip speed U^, the through -flow 
(axial) velocity U, and the work coefficient 8. The work coefficient is a measure of the 
loading on the fan. 
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In most high-Reynolds -number turbulent flows the characteristic frequency u>£ of 
the moving -frame turbulence correlation tensor is quite small compared with the fre- 
quency U c /£ of the turbulence observed in a fixed frame (Taylor’s hypothesis). In this 
case we can set equal to zero and the first term in equation (74) vanishes. Since 
the second term is only weakly dependent on it can be associated with noise gener- 
ated by the time unsteadiness of the turbulence. The second term therefore represents 
the generation of noise by the spatial nonuniformities and the convection of these non- 
uniformities by the mean flow. With u>£ equal to zero, equation (74) becomes 


pP 

m, n 


c nJ 3 U 
Crr ( 

-£ 2 ' 


| ^2 /2m?rA _ 2^p V 
V 6 A ) 


It 

2U 


co 


+ ai {i - x)u - 


27rplL 


A 


exp i 




r l w c 


L 


^ 2j7p ^r 
A p2 


U 


(78) 


The function W(Z/a) of Z exhibits a peak at Z = 0 whose width is equal to 2cr. Hence, 
it follows from equation (78) (with X - 1) and equation (75) that the spectrum will exhibit 
peaks at frequencies. 


U t U t 

co = 2irp — - = 2npB — - = pBft for p = 1, 2, . . . 

A 6 

where ft = 27rU t / 6 is the angular velocity of the rotor. Thus, equation (75) does indeed 
predict tone production at the multiple of the blade passing frequency. It can also be 
seen that the width of the tone is 4U c /Ir. The convection velocities of the eddies are 
probably nearly equal to the axial flow velocity U, which for most of the Quiet Fans is 
nearly equal to one -half the speed of sound. Typical fan spectra measured at Lewis 
exhibit tones which are of the order of 100 to 200 hertz in width. Hence, the turbulent 
eddies must be of the order of 1 to 2 feet in length in order to produce such tones. Since 
the most intense turbulence entering the fan is in the casing boundary layer, whose width 
is considerably less than this, it is fairly unlikely that the turbulence is producing the 
tones unless an extremely large amount of stretching of the eddies occurs. Thus, it has 
been postulated that when a nontranslating fan is operating out doors, large-scale turbu- 
lence shed from buildings and nearby structures is drawn into the fan and, due to the 
large, convergence of the mean streamline, the eddies become very elongated in the di- 
rection of flow. In any event the turbulence could still contribute to the broad sidebands 
which are invariably observed in the tones of fan spectra. 


29 



In figures 5 to 9, dimensionless acoustic spectra are plotted against dimensionless 
frequency for various values of the design a nd turb ulence parameters. Since the work 

coefficient 9 and the turbulence intensity ^u^/u 2 enter equation (75) simply as multi- 
plicative factors, we absorb these into the dimensionless acoustic spectrum. Since, as 
we have indicated, most of the subsonic Quiet Fans operated at the Lewis Research Cen- 
ter have an axial-flow Mach number of approximately 1/2, all the plots are drawn with 
M = 0. 45. It is also reasonable to assume that the convection velocity U c of the turbu- 
lent eddies is equal to the axial-flow velocity U. All figures are drawn for this case. 

In all plots except figure 9, co f is taken to be equal to 10‘ 4 Cq/A. This is essentially 
equivalent to putting it equal to zero in equation (74) (which amounts to using eq. (78)). 

The effect of varying the axial correlation length of the turbulent eddies is illus- 
trated in figure 5. It can be seen from figure 5(a) that when the correlation length Z in 
the transverse direction is as small as one -half the blade gap A and when the ratio r 
of the axial correlation length to the transverse correlation length is less than 1/2, there 
is virtually no evidence of tones. It is not at all unreasonable to suppose that much of 
the turbulence entering a full-size fan (most of which is in the casing boundary layer) is 
of this size. When either of the turbulence correlation lengths starts to become equal to 
the blade gap, the spectrum rapidly begins to form humps at the blade passing frequency 
and its harmonics. These humps are formed chiefly through a decrease in the sound 
power in the region lying between the harmonics of the blade passing frequency, with the 
level of the tones remaining relatively constant. The tones are quite broad until the 
correlation length is of the order of several blade gaps. Figure 5(b) shows that far nar- 
rower tones can indeed be produced by eddies which are 10 blade spacings long. 

It can also be seen from figure 5 that the power in the tone at the second harmonic 
of the blade passing frequency is nearly equal to (or in some cases greater than) the 
power in the fundamental. For most fans tested at Lewis the second harmonic is 
about 5 decibels lower than the fundamental. This is an indication that inlet turbulence 
may not be the chief producer of tones in these fans (at least through the quadrupole 
mechanism). 

The effect on the acoustic spectrum of varying the blade rotational Mach number M t 
is illustrated in figure 6. It can be seen from figure 6(b) that when the eddies are small, 
increasing M^. tends to increase the power radiated in the tones while leaving that ra- 
diated at intermediate frequencies unchanged. There also appears to be a tendency to 
shift the spectrum to higher frequencies as M t increases. 

In figure 7 the effect of varying the blade number B is shown, while figure 8 shows 
the effect of varying the blade height b. It can be seen from these figures that there is 
very little change in the shape of the spectra as these parameters vary and that doubling 
B or b roughly doubles (i. e. , increases by 3 dB) the sound power radiated at each fre- 
quency. But the cross-sectional area of the fan is proportional to Bxb. Hence, the 
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power radiated per unit cross-sectional area of the fan remains roughly unchanged as B 
or b is varied. This tendency is observed in the Quiet Fans tested at Lewis. 

Finally, figure 9 shows the effect of varying the moving-frame turbulence frequency 
cOj. In most cases which occur in practice, WjA/cq will be very nearly zero, and the 
value of 10" 1 used in figure 9 is much too large. However, this large value was chosen 
so that the effect of this parameter will show up. It can be seen that increasing 
merely fills in the spectrum between the tones. 


CONCLUDING REMARKS 

Inlet turbulence can be a dominant source of fan noise and at the high subsonic Mach 
number where current aircraft engine fans operate, it is quite possible that this noise is 
being generated through a quadrupole mechanism. But, for the high- subsonic -Mach - 
number regime, no quadrupole analysis of this noise source has been developed. The 
present report is written to fill this gap. The numerical results exhibit a number of 
features observed in actual fan spectra. However, they also show that, unless the tur- 
bulent eddies become extremely elongated as they enter the inlet of the fan, this mech- 
anism contributes mainly to the broadband noise spectrum. But there increasing evi- 
dence that such elongation does indeed occur when the fan is not translating and that as 
a result inlet turbulence is the dominant source of the pure tones. The present results 
may be helpful in resolving this question. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 7, 1974, 

501-04. 
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APPENDIX A 


EVALUATION OF 
For < 0 equation (68) becomes 




(? 2 + u t T > 


+ exp<- 



/MMA 

+ U t r) + — i (UcT+d) 

\ rj 


+ 




X 



+ U t r) + 



+ exp<- 



+ 


u tT ) 





Upon carrying out the integrations this becomes 




Similarly for ^ > 0 
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Hence, upon simplifying the results, we find that 
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APPENDIX B 


SYMBOLS 


^m, n 
A„ 


m, n 


0 




V H r H I 

l u> 

9m 

i 

k 

*Sn, n 

L i> 

l 

M 

M t 

< n 
m, n 

n 

P(cu) 

P ± (w) 
m, it ' 

&\u)) 

P 


defined by eq. (76) 

cross-sectional area of = Constant plane 

number of blades 

blade height 

defined by eq. (18) 

speed of sound 

(0 r + iMM 2 )//3 2 

viscous stress tensor 

function of t or turbulence generating function 
Green r s function 

rotor velocity integrals, eqs. (35) and (36) 
axial component of acoustic intensity spectrum 
imaginary part 

wave number, w/c Q 

propagation constant, Vk 2- /3 2 4 n 

rotor velocity correlations, defined in eq. (51) 

turbulence correlation length in direction transverse to flow 

axial Mach number, U/Cq 

rotor tip Mach number, U^/cq 

defined by eq. (68) 

normal direction to wall 

sound pressure spectrum 

modes contributing to pressure spectrum 

power spectrum 

pressure measured above ambient pressure 



Po 

R 

R ijkZ J 

R u, k 

R i,jk’ R ij 


ljkl 

ik 

0 

ik 

Re 


m- 

st 


s 

T 

T. t 
13 

t 

U 

u c 

U r 

U, 


U 


U 1 

V 

V, 


W(z) 

Wj 

x,x. 


ambient pressure 

mean radius of annular duct 

turbulence correlation tensors 

total velocity correlation 

moving -frame turbulence correlation tensor 

spatial part of moving-frame turbulence correlation tensor 
real part 

ratio of correlation length in axial direction to correlation length in trans- 
verse direction 

duct surface 

large time interval 

Lighthiirs stress tensor relative to moving frame 
time (observer) 

axial velocity of flow entering fan 
eddy convection velocity 
relative velocity, V u t 2 + u2 
rotor blade-row velocity 
turbulence velocity 

mean square axial turbulence velocity 
spectrum of axial acoustic velocity, v| 
rotor velocity components relative to blade 
fluid velocity, acoustic velocity, i = 1, 2, 3 
constant reference velocity (complex) 

v i - 6 u u 
-z 2 

e erfc(-iz) 

phase-locked rotor velocity field 
coordinate system fixed to blades 
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x { ,x 

Y 


y>y t 

z 

z 


Of 


a 


a 


~2 

a 


r^j 



0 

^r 

r m, n 

r P 
m, n 


r 


o 


r 

A 


a 

6 ij 

a(x) 

z v Z 
7, % ? 

0 


K m, n 


\n, n 
v(t) 


coordinates of observation point 

iJ)l.rj 2 ,(i?3 +U 3 )/2} 

coordinates of source point 
X 1 + i/3 r X 2 

t ? 1 D + ii?2 

same as 

defined by eq. (17), i * 1, 2, 3 
defined by eq. (58) 

defined by eq. (70) 
defined by eq. (69), i = 1, 2, 3 

^1 - M 2 

V 1 - (v c o> 2 

inner product of duct eigenfunctions 

defined by eq. (74) 

circulation 
stagger angle 
interblade spacing 
defined by eq. (40) 
transverse duct length, BA 
Kronecker delta 
delta function 

turbulence moving -frame coordinates 
rotor-locked coordinate system 

work coefficient of rotor, -Aw 2 /U^ 

eigenfunction for transverse modes defined in eq. (9) 

parameter which approaches unity 

defined by eq. (22) 

region of duct external to blades 



£. separation coordinate r/' - rj 
p 1 acoustic density fluctuation 

Pq average density 

r time (source) 

*m, n duct eigenfunctions 

q shaft rotational frequency, 2u\J ^/ 5 

°m, n defined ^ «!■ (24) 

n defined by eq. (57) 

w frequency 

characteristic frequency of turbulence 
Subscripts: 

i integer 1, 2, 3 

] integer 1, 2, 3 

k integer 1, 2, 3 

l integer 1, 2, 3 

m integer 0, ±1, ±2, . . . 

n integer 0, 1 , 2, . . . 

p integer ±1, ±2, . . . 

q integer ±1, ±2, . . . 

Superscripts: 

( ) complex conjugate or time average 

t measured relative to a coordinate system moving with flow 
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Dimensionless sound power, 10 log 10 









Figure 9. - Effect of turbulence frequency (moving framef on acoustic 
spectra. Number of blades, B, 30; ratio of blade height to interblade 
spacing, W A, 4; axial Mach number, M, 0.45; rotor Mach num- 
ber, M», 0.8; turbulence correlation length in direction transverse 
to flow, [, 1; eddy correlation length ratio, r, 1. 
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